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Abstract 

The hyperfine spin splittings in heavy quarkonia are studied using the re- 
cently developed renormalization group improved spin-spin potential which is 
independent of the scale parameter /i. The calculated energy difference be- 
tween the J/ijj and the r\ c fits the experimental data well, while the predicted 
energy difference AM p between the center of the gravity of l 3 -Po,i,2 states and 
the l l P\ state of charmonium has the correct sign but is somewhat larger 
than the experimental data. This is not surprising since there are several 
other contributions to AM p , which we discuss, that are of comparable size 
(~ 1 MeV) that should be included, before precise agreement with the data 
can be expected. The mass differences of the ip' — rj' c , T(IS') — r]b, T(25) — rj' b , 
and B* — B c are also predicted. 



PACS numbers: 12.39.Pn, 12.38. Cy, 12.38.Lg, 12.39.Hg 







I. INTRODUCTION 



The hyperfine splittings in heavy quark-antiquark systems can provide information about 
strong interactions or Quantum Chromodynamics (QCD) at low energies. Since the motion 
of the heavy quark and antiquark in heavy quarkonia is nonrelativistic, their dynamics can 
be well described by nonrelativistic potential models. The hyperfine splittings arise from 
higher order relativistic corrections and can be calculated using perturbation theory, given 
the appropriate spin-dependent potential. Recently, significant progress has been made in 
the theoretical study of the spin dependent potential M. In this work the spin- dependent 
potentials were derived from QCD first principles using the Heavy Quark Effective Theory 
(HQET) P|- The spin-dependent potential was separated into short distance parts involving 
Wilson coefficients and long distance parts which were expressed in terms of gauge invariant 
correlation functions of the color-electric and color-magnetic fields weighted by the Wilson 
loop path integral ||. If the tree level values for the Wilson coefficients are used the potential 
reduces to Eichten's and Feinberg's result 0. And using the one-loop values of the Wilson 
coefficients, also calculating the correlation functions to one-loop in perturbation theory, 
the spin-dependent potential at the one-loop level in perturbation QCD ||]|| is recovered. 
However, the leading logarithmic terms appearing in perturbative calculations were also 
summed up in Ref. [J]] using the Renormalization Group Equation (RGE) to obtain a scale 
independent result. Therefore, the spin-dependent heavy quark-antiquark potential derived 
in Ref. [p]] is scale-independent and thus improves upon and generalizes both Eichten's and 
Feinberg's result || and the one- loop perturbative result ||||. In addition, this improved 
result Jl| satisfies all the general relations among the different parts of the spin-dependent 
potential ||. In the following we use this improved, more general potential [|TJ to calculate 
the hyperfine splittings in the cc, cb, and bb systems. Specifically, we calculate the energy 
difference between the 3 Si and the 1 So states and the difference AM P between the center of 
the gravity of the 3 -Po,i,2 states and the 1 P\ state. 

First we note that the P-wave hyperfine splitting AM P in charmonium has been exper- 
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imentally determined to be —0.9 ± 0.2 MeV 0, which is not only much smaller than the 
splittings caused by the spin-orbit and the tensor interactions, but also the S'-wave hyper- 
fine splittings, which typically are 50 — 100 MeV. Naively, one might estimate the hyperfine 
splitting to be smaller than the spin-orbit and tensor splittings by the order of v 2 , where v is 
the relative quark- ant iquark velocity, or about 1/10 in charmonium. This interesting point 
has been studied previously §,1,0,0,0,0]- According to the Fermi-Breit formula, which 
follows from lowest order perturbation theory, the hyperfine spin splitting is proportional to 
the wavefunction at the origin, which vanishes for P-waves. However, one-loop corrections 
give logarithmic terms that are nonlocal and allow a non-zero contribution to the P-wave 
hyperfine splittings. Several previous calculations |]8|j9] JT0|JTT|| of AM p used only the one-loop 
perturbative spin-spin potential |4]|| and the results are remarkably close to the experimen- 
tal value. This agreement with the experimental value of AM P , taking into account only the 
one-loop contribution is surprising, since there are other contributions to AM P of similar 
size; for example, nonlocal terms coming from higher orders. We will discuss such effects 
below, as it is instructive to see how these higher order contributions could affect the results. 

In the following we will use the general formula for the spin-spin part of the 
renormalization-group-improved spin-dependent potential that was derived in Ref. |]J to 
calculate the hyperfine spin splittings in the cc, bb, and cb systems. Since the spin-spin 
potential is a short distance feature, perturbation theory can reliably be used in the cal- 
culation. Our result for the l 3 Si — l 1 ^ splitting between the J/ip and the 7] c agrees well 
with the experimental value and our predictions for the mass differences if)' — rj' c , T(IS') —f]b, 
T(2S) — T]' b , and B* — B c are reasonable. However, the contribution to the P-wave energy 
difference, AM p , between the center of the gravity of l 3 Po,i,2 states and the l x Pi state, while 
having the correct sign, is somewhat larger than the experimental data. That is, when the 
contributions of the leading logarithmic terms are summed up and included, the agreement 
with that data is not as good as when only the one-loop perturbative spin-spin potential, 
in which the leading logarithmic contributions are not summed up and included is used. 
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We will discuss the implications of these results in greater detail below and point out that 
there are several other contributions to the rather small energy difference AM p which esti- 
mates indicate are of the same order of magnitude as the spin-spin contribution. It therefore 
appears that the agreement of the one-loop perturbative result with the data is probably 
fortuitous. 

The following Sec. II is devoted to the calculational methods. In Sec. Ill we present our 
numerical results and in Sec. VI we discuss these results and our conclusions. 



II. CALCULATIONAL METHODS 

To calculate the hyperfme splittings in the heavy quark-antiquark systems we will use the 
spin-spin part of the renormalization-group improved general formula for the spin- dependent 
potential |J derived in the framework of HQET || . In the derivation the renormalized two- 
particle effective Lagrangian was first calculated to order 1/m 2 . Then, treating the terms 
of higher order in 1/m in the effective Lagrangian as perturbations, the four point Green's 
function on the Wilson loop [15| with the time interval T was calculated in the limit where 



m — ► oo first followed by T — > oo |16||J. In this limit, using standard perturbative methods, 
the large T behavior of the Green's function is of the form 

I oc e~ T < m > r \ (1) 

From Eq. (p]) e(m,r), the potential energy between the quark and the antiquark, can be 
extracted. Expanding e(m,r) in powers of 1/m each of the spin-dependent potentials can 
be factorized into a short distance part, involving Wilson coefficients, and a long distance 
part, which can be expressed in terms of correlation functions of the color-electric and color- 
magnetic fields weighted by the Wilson-loop integral. Using the notation of Ref. [I], the 
resulting spin-spin potential is 

S • S 

AH ss (m 1 ,m 2 ,r) = — ^ c 3 (/x, mi)c 3 (/x, m 2 )V i (fj 1 , r) - 6N c g 2 (n)d(fx)S 3 (r) , (2) 

6m\m2 1 J 



3 



where mi, 777.2, and Si, S2 are the masses and the spins of the heavy quark and the antiquark, 
respectively, \i is the renormalization subtraction point, iV c is the number of colors, and g s (p) 
is the running coupling constant. The Wilson coefficients 03(11,7x1) and d([x) were calculated 

respectively, and are 



in leading logarithmic approximation in Ref. |T3| and Ref. 

f <*s(lJ>) 



9_ 

2r, 



c 3 (/i,m) 



v a s (mi) 



(3) 



and 



N 2 -1 

= oat9. c 3 (m 2 ,mi)[l - Cg(/i,m 2 )] 



8iV| 
iV 2 - f /a 



,(m 2 ) \ 25 



f - 



a 



18 

25 



(4) 



8iv| V a s( m i)y V a «( m 2)/ 

In Eq. (|2]) V^/i, r) is the color magnetic- magnetic correlation function which can be expressed 
as 



dz / ^ / ^(S i (xi,^) J B i (x 2 ,y))/(l), 

-T/2 J—T/2 1 



(5) 



where (• • •) is defined by 



(•••} = J [dA^Tr jp exp (ig t) dz, 



exp(iS Y M(A)), 



(6) 



Here C(r, T) represents the Wilson loop [Tgl , P denotes the path ordering, and r = |xi — x 2 |. 

We emphasize that this is a general result for the hyperfine part of the spin- dependent 
potential to order 1/m 2 . It absorbs the short distance contributions to the potential into 
the coefficients 03(11,771) and d(fi) while the long distance contributions to the potential are 
contained in the correlation function V^(fi, r). Moreover, the result is independent of the 
factorization scale since the //-dependence in the coefficients cancels the //-dependence in 
the correlation function. The first term in the bracket in Eq. (0) is a nonlocal term while 
the second term is a local one which is generated by mixing with the first (nonlocal) term 
under renormalization. We note that if the coefficients are evaluated at tree level; i.e., 
c 3 (/i, m) = 1 and cf(//) = 0, the potential reduces to the Eichten-Feinberg result 0. And if 
these coefficients are expanded to order a s (fi) and the correlation function is also evaluated 



only to one-loop, the logarithmic terms in Eq. (|2]) then reduce to the one-loop spin-spin 
potential Therefore, this renormalization-group improved potential, Eq. (2), extends 

both Eichten's and Feinberg's result || as well as the one-loop perturbative potential |§||, 
containing each of these results as special cases. 

In a heavy quarkonium state the typical momentum transfer is of order mv , where v is 
the relative velocity of the heavy quark and the antiquark, and the typical size is of order 
\/{mv). In such a low momentum region the correlation function V±(ll, r) could in principle 
have nonperturbative contributions and should, therefore, be calculated using nonpertur- 
bative methods. However, since confinement in QCD is color-electrical, it is reasonable to 
expect the color-magnetic field to be predominately a short distance effect. Thus the color 
magnetic-magnetic correlation function V^li, r) should fall off quite fast when the distance 



r becomes large. This is confirmed both by lattice calculations |18j and by the experimental 
fact that AM p is empirically very small. If V^/i, r) had a significant long distance component 
AM p would be considerably larger, contrary to the data. We therefore can safely assume 
that the potential V^/i, r) is a short distance effect which can be calculated using pertur- 
bative QCD. Its perturbative expression can be obtained from the following arguments: As 
mentioned above, the result for AH ss (r) in Eq. (^) is /i-independent since the //-dependence 
in the coefficients cancels the //-dependence in the correlation function. However, to explic- 
itly demonstrate this cancellation to all orders one must calculate the correlation function 
to all orders, which is impossible to do directly. Fortunately, using the RGE, this can be 
done in the leading logarithmic approximation. In momentum space the Fourier transfor- 
mation of V4(/t,r), denoted by V4(/t,g), is dimensionless and is only a function of the two 
variables /i and q. It must, therefore, be a function of ln(g 2 //t 2 ) and these logarithms can 
be summarized using the RGE in the effective theory. Alternately, there is another simple 

approach: If we choose li — q all of these logarithmic terms vanish and only the tree level 
N 2 -1 

term — ^ — 9s(p)\i*=<i remains in V^/i, q) in the leading logarithmic approximation. Then 
all the nonlocal logarithmic terms having been absorbed into the coefficients c^(q, m) and 
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d(q) . Consequently, we find the hyperfine part of the spin-dependent potential in momentum 
space in the leading logarithmic approximation to be 

iV 2 -l 



AH ss (mi,m 2 ,q) = ^ 2 g 2 s (q) 



-c 3 (q } mi)c 3 (q } m 2 ) - QN c d(q) 



(7) 



where c 3 (q,m) and d(q) are given by Eqs. (||) and @. This final formula, Eq. (|7|), for the 
hyperfine spin-dependent potential, which we will use in our calculations of the hyperfine 
spin-splittings in the cc, cb, and bb systems, improves upon the one-loop perturbative cal- 
culation in two important respects: (i) it is independent of fi and (u) it includes the higher 
order logarithmic terms. 

To first order perturbation theory in AH SS the energy shift caused by AH ss (r) is 

AE = J dh^l^AHUr^t,^) (8) 

where \I/y z (r) is the nonrelativistic wavefunction of the bound state with total angular 
momentum / and z-component l z . For simplicity we suppress spin and color indices and 
retain only the space-dependent indices. Separating the radial part, u(r), we write \l/^ z (r) 

as 

^(r)= U (r)Y^(M), (9) 

where Y^ m (8,4>) are the standard spherical harmonics. Rotational invariance implies that 
AE is independent of l z . Averaging over l z and using properties of the spherical harmonics, 
AE can be expressed as 

AE = J ^-\u{r)\ 2 AH ss {r). (10) 
Taking the Fourier transform, in momentum space AE is given by 

AE = I (|^) A ^)' ( U ) 



where 



f d^T 1 r 

£(q) = J — e iq ' r |w(r)| 2 = - J drr sin qr\u{r)\ 2 . (12) 



Finally, doing the angular integration we have 

AE = i/ (13) 

which we will use to numerically calculate the hyperfine splittings. Of course, the radial 
wave function u(r) must first be obtained by numerically solving the Schrodinger equation 
with some chosen potential and then [Eq. (|12|)] can easily be calculated by using the 
Fast Fourier Transformation program. 



III. NUMERICAL RESULTS 



Using Eqs. (0), ([12]), and ( |i~3"D the hyperfine spin splittings for both the 5"- wave and the 
P-wave states were numerically calculated. The radial wavefunction was obtained by nu- 
merically solving the Schrodinger equation. For comparison, we used three popular potential 



models. One was the Cornell model |19| in which the potential has the form 



V(r) = -- + - 2 ., 
r a z 



(14) 



with 



m c = 1.84 GeV , m b = 5.18 GeV 
k = 0.52 GeV , a = 2.34 GeV 



(15) 



The second one was the logarithmic potential |2(| given by 



V{r) = -0.6635 GeV + (0.733 GeV) log(r • lGeV), 



(16) 



with 



m c = 1.5 GeV , m b 



4.906 GeV. 



(17) 



The third one was the improved QCD-motivated potential M with the form 



V{r) = - 



16vr 1 
25 rf(r) 



1 + 



2 7£ + § 462 In /(r) 



f(r) 



625 f(r) 



(18) 
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where f(r) was given by 



/(r)=21n 



' A ms + ( A n - a ms) ex P - t 15 (0.75A n - A m 



with 



A H A MS r 



m c = 1.478 GeV , m b = 4.878 GeV, 
An = 0.72 GeV, a = 2.59 GeV" 1 , 
C = 4.62. 



C 



(19) 



(20) 



To proceed with the calculation we also required an expression for the running coupling 
constant a s (q). The familiar RGE, one-loop result is 

47T 



a s (q) 



2 ' 



(21) 



A^ 

MS 



where bo = 11N C — 2Nf and Nj is the number of quark flavors. It is clear from Eq. fl2~l~|) 
that a s (q) contains a Landau singularity in the nonperturbative region when q 2 = A 2 ^ 
and becomes negative for q 2 < A. 2 -^. To avoid the resulting numerical ambiguities we first 
moved this singularity to q 2 = and used a modified form of a s (q) in the actual numerical 
calculations; namely, 

47T 



a s (q) 



b In 



(22) 



A^ 

< MS 



+ 1 



In the next section we shall discuss alternative approaches, the sensitivity of the results, and 
their implications. The value of Aj^g was taken to be 200 MeV and 250 MeV in the numerical 
calculations, which is within the experimental range A^ = 195 + 65 — 50 MeV |21| . Our 
numerical results for the three potentials are presented in Tables I, II, III, respectively. For 
comparison we have also included the results for the 2S and 2P states. The main features 
of these results can be summarized as follows: 



The results are //-independent, as they must be. 
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The calculated energy difference between the J/^> and the r] c mesons is quite close to 
the experimental value for all three potentials. 

For each of these three potentials we predict the energy difference between \J r ' and rj' c 
to lie within the range 55 — 80 MeV. 

For the bb system there are significant discrepancies between the Cornell model, with 
the parameters given by Eq. QI5D , and the other two models for the S-states. Since 
the Cornell model, with these parameters, does not predict the bb spectrum very well, 
the results calculated in the other two models are probably better predictions for the 
energy difference between the T(1S) and the rjb (35 — 50 MeV) and between the T(2S) 
and the r]' b (20 MeV). 

The predicted energy difference between B* and B c meson is in the range 40 — 70 MeV 



from all three of these models, which is consistent with previous results [13 



The calculated value of AM p = E(l 3 Pj) — E(l l Pi) for the charmonium IP states is 
in the range of —4 to —6 MeV, which has the same sign but is several times larger 
than the experimental value of —0.9 ± 0.2 MeV 0. This is not surprising since there 
are several other contributions to AM P which estimates indicate are comparable in 
magnitude to the contribution coming from the hyperfine spin-spin interaction, H ss . 
In fact, it is surprising that the prediction from only the one- loop spin-dependent 
potential is quite close to the experimental data. We discuss these other contributions 
in the next section. 



IV. DISCUSSION AND CONCLUSIONS 

We have calculated the hyperfine spin splittings in the cc, bb, and be system using the 
RGE improved perturbative spin-spin potential The results for the hyperfine splittings 
of the S"-wave states agree with the J — r] c measured splitting and the prediction for 
splitting T — r]b is reasonable. However, the contribution to AM P = E( 3 Pj) — E{}Pi) for 
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the charmonium P-wave states is somewhat larger than the experimental data , although 
it agrees in sign. That is, after summing up the leading logarithmic terms and including 
them in the perturbation calculations, the agreement with the data is not as good as the 



<M?) = ^ jj = a s {fi) 

1 - 73a s (/i) ln- 2 



-" 2 > (23) 



m=l 



47T q 2 



one-loop calculations [p|,p|, [T0] , |TT]| . In order to illustrate this clearly, we can expand a s (q) in 
terms of and truncate it at some finite order. In our final formula, Eq. (|?p, we used 

the expansion for a s (q), 

— a s (// In — 

47r q z 

and truncated at several choices of n. Specifically, we repeated the numerical calculations 
for the improved QCD motivated potential |§ for n — 1, 2 and 4, choosing the scale /i, 
now to be \i = 1.5 GeV, 4.0 GeV, and 2.5 GeV for the cc, bb, be systems, respectively. The 
numerical results are presented in Tables IV, V, and VI corresponding to n = 1, 2, and 
4, respectively. For comparison, we also presented the results obtained using the complete 
one-loop hyperfine potential in Table VII. Comparing Table IV and Table VII we see 
that AMp for n — 1 is quite close to the complete one-loop result. However, from Table V 
and Table VI we see that the predicted values of AM p are about 60%-80% and 150%-200% 
larger than when terms up to order 2 and order 4 are kept in the expansion of a s (q), Eq. (P3[). 
We note that we also repeated these calculations for the logarithmic potential f2(J and the 
Cornell potential ||19|| . All three potentials predicted similar values for AM p . This clearly 
indicates that the nonlocal logarithmic terms from high loop perturbative calculations are 
quite important. In fact, even using the RGE to sum up these logarithmic terms does not 
allow one to understand the experimental value of AM p , indicating that the success of the 
one-loop calculations |]8|j9| J10|JTT[1 was probably fortuitous. In fact, there are several additional 
contributions that are possibly comparable in magnitude. These include the following: 

The contributions of the spin-orbit and and tensor potentials in the second order of 
perturbation theory: These contributions to AM P only cancel to first order in perturbation 
theory. However, according to the power counting rules introduced in Ref. [E!3], the spin- 



orbit and tensor potential potentials shift the energies of the P-wave states by an amount 
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of order mv 4 in first order, which indeed cancel in AM P , but they do make a contribution 
to AM p of order mv 6 in the second order of perturbation theory. This estimate is several 
MeV for the P-wave charmonium states, and therefore should not be neglected. 

Higher dimensional operators: Unlike the dimension-six operators, these give non-zero 
contributions to AM p even at tree level. Compared to the one-loop contribution, these are 
suppressed by v 2 but enhanced by aj 1 and v 2 /a s ~ 1 in charmonium. 

The color-octet S'-wave component in P-wave quarkonia states [f^]: This component of 
the wavefunction receives a tree- level contribution from the local term <5 3 (r) in the spin-spin 
potential. This contribution too could be of order v 2 /a s ~ 1 compared to what has been 
calculated. 

Next-to-leading order perturbative contributions from the two-loop potential: These are 
suppressed by order a s , but since a s is not a very small quantity in charmonium, one cannot 
dismiss the possibility that this contribution could be significant. 

Before comparing with the experimental value of AM p in charmonium, which is only 
about 1 MeV, all the above contribution should be included since they are possibly compa- 
rable in magnitude. In the bb case these effects are less important and one can expect the 
perturbative calculations the bb system to be more reliable, although AM p is smaller, also. 

Finally, to explore the sensitivity of our results to the location of the Landau singularity 
in a s (q) we replaced the expression, Eq. (20), by 

a s (q) = 7—2 s- • (24) 

6oln xt + A2 

V MS / 

and varied A 2 . The results for the S- wave hyperfine splitting were not sensitive to A 2 and 
only for large A 2 did AM p significantly decrease. To fit AM p to the measured value required 
A 2 quite large, about 16, clearly out of the perturbative region. 

We thank Eric Braaten, Yu-Ping Kuang for their instructive discussions and comments. 
This work was supported in part by U.S. Department of Energy, Division of High Energy 
Physics, under Grant No. DE-FG02-91-ER4086. 
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Table I. The hyperfine spin splittings in MeV predicted by Eq (|7|) with Cornell 
potential |19| 

cc bb be 



Ams (MeV) 


200 


250 


200 


250 


200 


250 




117.1 


128.3 


97.7 


104.3 


67.0 


71.4 


E(2 3 S 1 ) - E(2 1 S ) 


75.3 


82.5 


39.6 


42.2 


37.7 


40.3 


E^P^-E^P,) 


-4.8 


-7.0 


-3.0 


-4.1 


-4.0 


-5.6 


E(2 3 Pj) - E{2 1 P l ) 


-3.6 


-5.2 


-2.1 


-2.9 


-3.0 


-4.2 



Table II. The hyperfine spin splittings in MeV predicted by Eq (|7|) with Logarithmic 
potential |20| 

cc bb be 



Ams (MeV) 


200 


250 


200 


250 


200 


250 


E(1 3 S 1 )-E(1 1 S ) 


106.1 


117.1 


36.0 


38.2 


40.6 


42.6 


E(2 3 S 1 ) - E(2 1 S ) 


54.2 


59.7 


18.7 


19.8 


21.2 


22.3 


E^P^-E^P,) 


-5.4 


-7.8 


-3.4 


-4.5 


-4.5 


-6.4 


E(2 3 Pj) - E{2 1 P l ) 


-2.6 


-3.8 


-2.1 


-2.8 


-2.8 


-4.0 



Table III. The hyperfine spin splittings in MeV predicted by Eq (0) with Improved- 
QCD motivated potential [§] 

cc bb be 



Ams (MeV) 


200 


250 


200 


250 


200 


250 


^l 3 ^)-^! 1 ^) 


107.9 


119.1 


44.6 


47.6 


43.4 


45.7 


E(2 3 S 1 ) - E(2 1 S ) 


68.5 


75.6 


20.9 


22.4 


25.2 


26.7 


E^P^-E^P,) 


-4.6 


-6.7 


-2.7 


-3.7 


-3.7 


-5.3 


E(2 3 Pj) - E{2 1 P l ) 


-3.4 


-5.0 


-1.8 


-2.5 


-2.6 


-3.8 



14 



Table IV. The hyperfine spin splittings in MeV predicted by Eq (^) with n = 1 for 
the Improved-QCD motivated potential || 

cc bb be 



Ams (MeV) 


200 


250 


200 


250 


200 


250 




111.8 


125.4 


45.7 


49.2 


45.2 


48.6 


E(2 3 S 1 ) - E(2 1 S ) 


71.3 


80.1 


21.5 


23.2 


26.3 


28.4 


E^P^-E^P,) 


-1.5 


-1.9 


-1.0 


-1.2 


-1.2 


-1.5 


E(2 3 Pj) - E{2 1 P l ) 


-1.0 


-1.3 


-0.7 


-0.8 


-0.9 


-1.2 



Table V. The hyperfine spin splittings in MeV predicted by Eq (^) with n = 2 for 
the Improved-QCD motivated potential |§ 

cc bb be 



Ams (MeV) 


200 


250 


200 


250 


200 


250 


E(1 3 S 1 )-E(1 1 S ) 


110.5 


123.4 


45.2 


48.5 


44.7 


47.8 


E(2 3 S 1 ) - E(2 1 S ) 


70.6 


79.1 


21.3 


23.0 


26.1 


28.2 


E^P^-E^P,) 


-2.5 


-3.4 


-1.6 


-2.0 


-2.0 


-2.7 


E(2 3 Pj) - E{2 1 P l ) 


-1.6 


-2.2 


-1.0 


-1.3 


-1.4 


-1.8 



Table VI. The hyperfine spin splittings in MeV predicted by Eq (|7|) with n = 4 for 
the Improved-QCD motivated potential |§ 

cc bb be 



Ams (MeV) 


200 


250 


200 


250 


200 


250 


^l 3 ^)-^! 1 ^) 


108.8 


120.4 


44.8 


47.9 


43.9 


46.5 


E(2 3 S 1 ) - E(2 1 S ) 


69.3 


76.8 


21.0 


22.7 


25.6 


27.4 


E^P^-E^P,) 


-4.0 


-6.2 


-2.3 


-3.0 


-3.1 


-4.3 


E(2 3 Pj) - E{2 1 P l ) 


-2.8 


-4.3 


-1.4 


-1.9 


-2.1 


-2.9 
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Table VII. The hyperfine spin splittings in MeV predicted by the complete one-loop 
spin-spin potential [|]|| with Improved-QCD motivated potential || 

cc bb be 



Ams (MeV) 


200 


250 


200 


250 


200 


250 


E(1 3 S 1 )-E(1^S ) 


127.7 


145.7 


50.1 


54.3 


39.2 


41.5 


E(2 3 S 1 ) - E{2 1 S ) 


81.5 


93.3 


23.4 


25.6 


22.6 


23.9 


E^P^-E^P,) 


-1.2 


-1.5 


-0.6 


-0.7 


-0.6 


-0.7 


E(2 3 Pj) - E(2 1 P 1 ) 


-1.1 


-1.3 


-0.5 


-0.6 


-0.5 


-0.6 
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